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Abstract. Dispersive approaches provide model independent description of meson-meson scatter- 
ing. We first review here the use of dispersion relations to obtain a model independent unitarization 
of Chiral Perturbation Theory amplitudes, that establish the existence of light scalar mesons and 
whose leading 1 /N c behavior suggest they ahve a non qq dominant component. We also review the 
, forward dispersion relation checks on conflicting experimental data and the resulting very precise 

nn scattering scattering amplitudes. 
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(Nj ; Over the last years there has been a renewed interest in meson-meson scattering. 

The reasons are that its low energy description is relevant for understanding the QCD 
vacuum, and that the description below 1.2 GeV requires the existence of some poles 
in the amplitudes, related to light scalar mesons. Such resonances are the subject of a 
strong debate, first on their existence, that seems to be settling, but also on their nature. 
Data on meson-meson scattering is obtained indirectly from other processes involving 
qj \ nucleons or the decays of other heavy mesons. Thus, very often the existing data on 

meson-meson amplitudes are extracted with strong model dependent assumptions (one- 
>• . pion exchange, absorption models, rescattering, pole extrapolations to define the initial 

meson-meson state, etc.), and are therefore plagued with systematic uncertainties much 
larger than statistical errors. There are also K — > irnlv experiments, known as A74 
decays, where the pion phase shifts are extracted in a particular combination of isospin 
and 1, free of the previous systematic uncertainties and yield very precise and reliable 
results, but limited to invariant masses smaller than the kaon mass. 

Concerning the theoretical description, we have two model independent approaches. 
On the one hand, we have Chiral Perturbation Theory (ChPT), which is the effective 
Lagrangian of QCD, written as an expansion in masses and derivatives of pions, kaons 
and etas, which are the Goldstone bosons of the spontaneous chiral symmetry breaking 
of QCD. The only caveat to this systematic expansion is that it can only be applied 
at low energies. On the other hand, it is possible to use the usual S-matrix constraints 
of causality, analyticity, unitarity, crossing, etc... to write dispersion relations for the 
different meson-meson channels. We review here how we have recently applied this 
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approach to check the consistency of different data sets and to obtain a precise pion-pion 
scattering amplitudes in the whole energy range. Of course, both approaches can also 
be combined to obtain a model independent description of data that also incorporates 
the low energy chiral symmetry constraints (see. H. Leutwyler's talk on Roy equations 
in this conference). We review here how this has been done by means of the Inverse 
Amplitude Method (IAM), which has the advantage that only uses ChPT input in the 
dispersive integrals and therefore allows for a relation between the resulting fits and 
QCD. Since the IAM generates the light resonances that appear in scattering, it is then 
possible to study their nature in terms of QCD parameters, like the number of colors. 

THE INVERSE AMPLITUDE METHOD FROM DISPERSION 

THEORY 

We will review here how the one-channel Inverse Amplitude Method (IAM) [Q], 0, [H] 
for pion-pion scattering is obtained just by using ChPT up to a given order inside 
a dispersion relation. There are no further assumptions and therefore the approach is 
model independent and provides an elastic amplitude that satisfies unitarity and has the 
correct ChPT expansion up to that given order. 

To fix ideas, let us consider the ChPT series for a pion-pion scattering partial wave 

(2) (4) 

amplitude of definite isospin / and angular momentum J, namely, tjj = t u + f 7 y + ... 
where t 2 = 0(p 2 ), t$ = 0(p 4 ) and p stands for the pion mass or momentum. For the 
complete partial wave tjj(s), it is possible to write a dispersion relation 

2 , s 3 r lmt u (s')ds' 
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that, for convenience, we have subtracted three times. Note we have explicitly written 
the integral over the right hand cut (or physical cut, extending from threshold, s t h to 
infinity) but we have abbreviated by LC the equivalent expression for the left cut (from 
to — oo). We could do similarly with other cuts, if present, as in the TlK case. 

We can also write dispersion relations for f ( 2 ) and t^ 4 \ but remembering that is a 
pure tree level amplitude and it does not have imaginary part nor cuts: 

(2) (4) ? s 3 r ImJrV (s')ds' , (a) s 

tif=a + a lS , ti ) =b + b lS + b 2S 2 + - -^LL2_ +I£: ( f W). (2 ) 

We now recall that unitarity, for physical values of s in the elastic region implies: 
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lmt u = o\t u \ l => Im— = -o => t IJ = —— l — — , (3) 



where a = 2p/^fs. Therefore, the imaginary part of the inverse amplitude is exactly 
known in the elastic regime. We can then write a dispersion relation like that in (QQ) but 

now for the auxiliary function G = (tff) 2 /t/j, i.e., 

G( S ) = G + G lS + G 2S 2 + - r ^f is ' )dS ' +LC(G)+PC, (4) 

% Js th s'*{s' — s — ie) 



where now PC stands for possible pole contributions in G coming from zeros in t[j. 

(2) 

It is now straightforward to expand the subtraction constants and use that Im^ y = 
and Im?jy = o\tff\ 2 , so that ImG = — Imfjy . In addition, up to the given order, 
LC(G) ~ —LC(tjj^), whereas PC is of higher order and can be neglected. Thus 



We have thus arrived to the so-called Inverse Amplitude Method (IAM): 
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that provides an elastic amplitude that satisfies unitarity and has the correct low energy 
expansion of ChPT up to the order we have used. It is straightforward to extend it to other 
elastic channels or to higher orders [Q]. Note also that, by looking at ©, it seems that 

it can also be derived by replacing Refjy ^ by its 0(p 4 ) ChPT expansion. But, strictly 
speaking, © is only valid in the real axis, whereas our derivation allows us to consider 
the amplitude in the complex plane, and, in particular, look for poles of the associated 
resonances. Actually, already ten years ago |30, with the single channel IAM we were 
able to generate poles for the p(770), ^*(892) and most interestingly, the controversial 
o (also called /o(600)), without any model dependent assumptions. 

Note that the above one-channel IAM derivation is model independent, and that 
contrary to a wide belief in the community contains a left cut and respects crossing 
symmetry up to, of course, the order in the ChPT expansion that has been used. 

The confusion may come from the fact that the IAM has also been applied in a coupled 
channel formalism, for which there is still no dispersive derivation, and sometimes with 
further approximations. Indeed one can arrive to © in a matrix form, ensuring coupled 
channel unitarity, just by expanding the real part of the inverse T matrix. For the coupled 
channel case different approximations to ReT 1 have been used: 

• The fully renormalized one-loop ChPT calculation of ReT ' 1 provides the correct 
ChPT expansion in all channels, also with left cuts approximated to 0(p 4 ) [0,[5l]. Indeed, 
using ChPT parameters consistent with previous determinations within standard ChPT, 
it was possible J5|,@] to describe below 1.2 GeV all the scattering channels of two body 
states made of pions, kaons or etas. Simultaneously, this approach [60 generates poles 
associated to the p(770), AT* (892) vector mesons, and the /b(980), a (980), o and K 
(also called ^o(800)) scalar resonances. 

• Originally |7j], the coupled channel IAM was used neglecting the crossed loops 
and tadpoles. This approach is considerably simpler, and although it is true that the 
left cut is absent, its numerical influence was shown to be rather small, since the meson- 
meson data are nicely described with very reasonable chiral parameters and generates all 
the poles enumerated above. Let us remark that this approximation keeps the s-channel 
loops but also the tree level up to 0(p 4 ), and that this tree level encodes the effect of 
heavier resonances, like the rho. Thus, contrary to some common belief, this approach 
still incorporates, for instance, the low energy effects of t-channel rho exchange. 

• Finally, if one is interested in describing just the scalar meson-meson channels, 
it is possible to use just one cutoff (or even a dimensional regularization scale) that 



numerically mimics the combination of chiral parameters that appear in those scalar 
channels. This has become very popular, even beyond the meson-meson interaction 
realm, due to its great simplicity but remarkable success [8|]. 

NATURE OF LIGHT SCALARS FROM UNITARIZED CHPT 

One of the big advantages of the unitarization approaches described in the previous 
sections is that when they use the fully renormalized ChPT amplitudes, they therefore 
have the correct chiral and flavor symmetry structure, including both the spontaneous 
and explicit symmetry breaking. Furthermore we also have the correct dependence on 
QCD parameters like the number of colors, which is of particular interest, since there 
are sharp predictions on how the mass and width of qq resonances should behave in a 
large N c QCD expansion. In particular, M ~ 0(1) whereas T ~ 0(1 /N c ). 

The \/N c leading behavior of all ChPT parameters is known and model independent, 
so that they can be varied accordingly to study [9] the N c dependence of the amplitudes. 
In particular, we can study the N c behavior of all the poles generated in the IAM. It has 
been shown that both the mass and width of the vector mesons generated with the IAM 
follow remarkably well the expected qq behavior. However, light scalar s do not follow 
a dominant qq N c behavior [9] , at least for N c not too far from real life, N c = 3. In figure 
1 we illustrate these two different behaviors calculated at with the 0(p 4 ) IAM, for the 
p AT* (892), and for the o, and K, whose poles can be obtained with the one-channel 
IAM, and therefore in a model independent way. Similar plots can be found for the other 
scalars in BSD. 

At this point it is worth noting that two-meson loop diagrams are subdominant at 
large N c . Indeed, the above results imply that for vectors the meson loops play a small 
role in the cancellations that lead to poles in the amplitude, whereas for scalars such loop 
diagrams are very important at N c = 3. Since these diagrams become smaller and smaller 
one could wonder about the influence of higher order effects in scalars. Thus recently, 
we have performed ft 1 Oil the full two-loop IAM analysis of the o channel confirming 
that, as it happened for the 0(p 4 ) amplitude, close to N c = 3 the sigma behaves rather 
differently than expected for a qq state, but that as the loop diagrams are suppressed, a 
subdominant qq behavior is recovered at larger N c . Remarkably, this qq behavior arises 
slightly above 1 GeV, as it can be seen in Figure 1 (right plot). This seems to support 
the suggestion rfTTIl that there is a non-qq scalar nonet below roughly 1 GeV, and another 
qq nonet above, but using a mode independent framework based on dispersive integrals, 
ChPT and the large N c QCD behavior. 

PRECISE AMPLITUDES FROM DISPERSION RELATIONS AND 

ROY EQUATIONS 

In the previous approach we used ChPT inside the integrals of partial wave dispersion 
relations. This has the advantage that we can control all the parameters of our input, and 
allow for their variation in order to understand the nature of poles or the dependence on 
certain QCD parameters. However ChPT does not necessarily give a good description of 
the integrands at high energies. Thus, even though we performed several subtractions to 
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FIGURE 1. Leading N c behavior of resonance poles generated IAM. Left and center columns: with 
one-loop 0(p 4 ) ChPT Both the p and K* poles in the imaginary plane behave as expected for qq states, 
namely M ~ 0(1) and Y ~ 0{\ /N c ), respectively. In contrast, the light scalars a and K do not behave 
predominantly as qq. Right: with two-loop 0(p 6 ) ChPT. We plot the mass M and width F evolution with 
N c calculated from the pole position, normalized to the N c = 3 case. Once more the dominant behavior is 
not that of a qq but a subdominant qq behavior emerges at larger N c around 1 GeV in mass. 



suppress the high energy regime we cannot claim to have very precise amplitudes, but a 
qualitative and most likely a semi-quantitative description. 

There are however other approaches where one can use directly the data inside the 
integrals. Of course, now it is much harder, if possible at all, to interpret changes of 
parameters in terms of QCD, but we can get extremely precise results for the amplitude 
and other observables like scattering lengths, poles, etc... Recently, ChPT constraints 
and data have been included in single channel dispersion relations for %% and kK partial 



waves QX2|] , confirming the existence of poles for a and K resonances, introducing some 
cutoffs on the dispersive integrals. For partial waves, the left cut is always a very delicate 
issue, because the large t behavior is not well known. Two possible model independent 
approaches to this problem have been given in the literature: One is to rewrite the left 
cut integrals in terms of a coupled set of integral equations relating different partial 
waves, known as Roy equations for %% scattering and Roy-Steininger equations for %K 



scattering, which have been recently used to obtain precise determinations of the o IU3I1 
and K jl4 \ poles, respectively. 

Here I will comment on the other approach, namely, to use Forward Dispersion 
Relations (FDR), that is, to avoid using partial waves and use full amplitudes setting 
t = 0. In a recent analysis we have used the following set of dispersion relations that 
form a complete isospin set: by choosing either F = Fqq or F = Fq + in 

F(,)-F(4M 2 ) = ^- 4M - ) PP f ds> (^-4A^)ImiV) 



we have two dispersion relations which are twice subtracted . Thus, the weight of the 
high energy part is quite suppressed, indeed as much as in Roy equations, but with the 
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FIGURE 2. The data fits obtained in 1161 11511 satisfy remarkably well the complete set of FDR and 
Roy Equations below KK threshold and fairly well above. In this plots we compare the result of the 
direct parametrization (direct or "in" curves) versus the integral representation (dispersive or "out"). We 
emphasize these curves come from fits to data, without FDR or Roy equations as constraints, so that the 
agreement within uncertainties between continuous and dashed lines is even more remarkable. 



advantage of having always positive contributions to the integrand, a fact that makes 
them much more precise. In addition, by setting s = 2M%, and F = Fqq, we find two sum 
rules important to fix the Adler zeros. Finally, for the t-channel exchange of isospin 1, 
which does not require subtractions, we use, 

F ( - I '- 1 >(s,0) = £P.P. / ds'— — — v ' (8 

v ' % UmI (s'-s)(s' + s-4Ml) 

At threshold this is known as the Olsson sum rule. 



First of all, we have found 111511 that only a few sets of data available in the literature, 
satisfy reasonably well the above forward dispersion relations within errors. Some of 
the most widely used experimental phase shift determinations indeed fail, on average 
for each FDR, by more than 2 standard standard deviations. The sets that satisfy FDR 
better are those closer to the fit using just the low energy data on decays commented 
in the introduction. All these sets have a marked "hunchback" in the 400 to 900 MeV 
region that undoubtedly is caused by the o pole. We show in Figure 2 (upper row) 
the comparison between the FDR left hand sides ("direct" calculation) and FDR right 
hand sides ("dispersive" calculation). Note that the agreement below KK threshold 
is astonishing given the tiny uncertainties and that only data has been fitted, not the 
FDR themselves. We have recently checked [18] that our direct fits to data also satisfy 
remarkably well Roy equations, as shown in Figure 2 (lower row). 



In lfl5ll we went further than this and improved the kk amplitudes by constraining the 
fits to different partial waves to satisfy the FDR below 950 MeV. In this way we obtained 
an even more precise representation that lies not too far from the best fits to data only, 
precisely those showing the o "hunchback" in the SO wave. 

Compared with the solution from Roy equations in [fTvh . the phase shift solutions 
that satisfy better the FDR representation in Ill5n agrees within errors in the low energy 
region for the S and P waves, but the SO phase shifts above 400 MeV lie about 2 standard 
deviations higher than in Ill7ll . We also have a discrepancy of about two standard 
deviations in the D wave scattering lengths and the Regge residue of the p . 

At present, and using our very recent improved fits above KK threshold, we are 
working on an even more precise KK amplitude by constraining our fits to satisfy 
simultaneously FDR and Roy equations, together with Froissart-Gribov and other sum 
rules. This parametrization will be used to obtain a precise determination of the o pole 
as well as of threshold parameters to determine the low energy ChPT constants. 

In summary, I believe that, dispersive approaches combined with ChPT are the most 
powerful technique at our disposal to describe precisely and understand in terms of QCD 
the meson-meson interactions and light scalar spectroscopy. 
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